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We demonstrate that single-photon as well as biphoton revivals are possible in a new class of 
dynamic optical systems-the so-called Glauber-Fock oscillator lattices. In these arrays, both Bloch- 
like oscillations and dynamic delocalization can occur which can be described in closed form. The 
bunching and antibunching response of path-entangled photons can be pre-engineered in such cou- 
pled optical arrangements and the possibility of emulating Fermionic behavior in this family of 
lattices is also considered. We elucidate these effects via pertinent examples and we discuss the 
prospect of experimentally observing these quantum interactions. 
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The prospect of manipulating and engineering quan- 
tum states has become an issue of great importance 
within the framework of quantum information and com- 
putation [l|, Along these lines, several physical set- 
tings have been envisioned as viable avenues to achieve 
this goal. Among them, one may mention trapped-ion 
arrangements and optical lattices as well as spin systems 
and quantum dots [3|, While the list of such possi- 
bilities keeps increasing with time, quantum optics has 
so far provided a versatile platform where such ideas can 
be experimentally realized and tested. As previously in- 
dicated, in optics, quantum information processing can 
be achieved entirely linearly, using simple passive com- 
ponents like beam splitters and phase shifters along with 
standard photodetectors and single-photon sources [5]. In 
this same optical realm, quantum entanglement can arise 
as a natural byproduct of photon interactions-a clear 
manifestation of their particle- wave duality. Perhaps, 
nowhere this process is more apparent than in the so- 
called Hong-Ou-Mandel two-photon interference effect [6]. 
In this latter configuration, photon entanglement is made 
possible via quantum interference- afforded after scatter- 
ing from a beam-splitter. Lately, optical arrays of evanes- 
cently coupled waveguides have been suggested as a pos- 
sible route toward the implementation of multiport sys- 
tems with moldable quantum dynamics [7]. The flexibil- 
ity offered by such compact and often miniaturized op- 
tical N X N configurations is made possible by the ex- 
ceptional control achievable these days in microfabrica- 
tion techniques [s*, ^9^. In this regard, Bloch oscillations 
of NOON and W entangled states as well as quantum 
random walks have been theoretically considered and ob- 
served in such arravs [lQl - [l2j . In addition, the evolution of 
quantum correlations in both periodic and random (An- 
derson) lattices has also been investigated [13, 14]. The 
possibility of classically emulating Jaynes-Cummin^s sys- 
tems on such lattices has also recently been proposed |15|. 
The question naturally arises as to whether such multi- 
port array systems can be utilized as a means to manip- 



ulate quantum states of light. 

In this Letter we investigate the propagation dynamics 
of non-classical light in a new class of dynamic pho- 
tonic systems-the so-called Glauber-Fock oscillator lat- 
tices. We demonstrate that Bloch-like revivals and dy- 
namic delocalization effects can naturally occur in spite 
of the fact that the structure itself is semi-infinite and 
not periodic. Interestingly, these interactions can be 
described in closed form, from where one can analyti- 
cally deduce the turning points of these quantum oscilla- 
tions. More importantly, the bunching and antibunch- 
ing response of path-entangled biphotons can be pre- 
engineered in such coupled optical arrangements. Emu- 
lating Fermionic dynamics in such arrangements are also 
considered and compared to those expected from bosonic 
systems in these same arrays. Finally the possibility of 
experimentally realizing such Glauber-Fock oscillator lat- 
tices is discussed. 

We begin our analysis by considering a semi-infinite 
Glauber-Fock oscillator array consisting of evanescently 
coupled waveguides. In this arrangement the coupling 
coefficients among neighboring channels vary with the 
square root of the site index, i.e., Ck,k-\-i oc \/n + ijof. 
For generality, we also allow this coupling to depend 
on the propagation distance in this lattice, in which 
case C/c,/c+i (X f{Z)y/n + 1 where f{Z) is an arbitrary 
real function. In addition we also assume that the 
propagation constant (local eigenvalue) of each waveg- 
uide element varies linearly with the site position. In 
essence, in this arrangement the refractive index is lin- 
early increasing-in a way analogous to that of an exter- 
nally biased crystal in solid state physics. Starting from 
these premises, one can show that in this class of arrays, 
the Heisenberg equation of motion for the creation op- 
erator of a single photon in waveguide mode k is given 
by: 

i^-f{Z)[Vk+lai^, + Vkai_,)-\kai = (1) 
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In the above equation, Z represents a normalized prop- 
agation distance, and A is a real constant associated with 
the strength of the aforementioned linear index change 
among adjacent sites. We emphasize that unlike stan- 
dard infinite Bloch oscillator arrays [16rJ^], the proposed 
structure is semi-infinite and asymmetric, e.g. the waveg- 
uides are no longer equidistant. As we will see later, these 
additional degrees of freedom may enable one to observe 
Bloch-like oscillations even in the neighborhood of the 
array boundary {k = 0). This is in contradistinction to 
the well-known Dunlap-Kenkre system, the only other 
integrable oscillator lattice [16]. 

In general, the quantum dynamics in this Glauber- 
Fock oscillator array can be described through the evo- 
lution matrix T relating the input-output states, i.e.. 



4(0) = £T,*„(Z)at(Z). 



(2) 
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In Eq.(2), T^^ (Z) represents the Hermitian conjugate 

of the (/c, n) element of T matrix or unitary transforma- 
tion. We would like to emphasize that in the present 
case, the evolution matrix cannot be simply obtained 
from ex.p{—iZH) since the Hamiltonian of the problem 
is Z (or time) dependent [19]. Yet, in spite of this com- 
plexity, one can show that the evolution elements ^ of 
the system can be obtained in closed form. Starting from 
Eq.(l), one can show that these elements are given by: 
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and A (Z) = J,^ [jf e^A(z'-z'')/ (^/) dZ') f {Z") dZ" 
B{Z) = -ij^e-'^^'f{Z')dZ', and C {Z) = 



-iXZ 



B*{Z). In Eqs.(|3]) (^) represents the 
associated Laguerre polynomials. In addition, Eqs.(|3]) 
imply that Xl^o l^fc,n|^ = 1, in agreement with the fact 
that T is itself a unitary transformation. In order to 
gain insight into the quantum dynamics in this class of 
arrays, let us first consider the case where only a single 
photon is launched into the kih waveguide element. We 
begin by analyzing here the simplest possible scenario 
where the Hamiltonian of the system is Z— independent, 
thus / (Z) = 7, where 7 is a real constant. Under these 
conditions, the evolution matrix elements are given by 



Tk n 



n < k 
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where 6 = ij^X'^ + 7^ [exp(-zAZ) - 1] A-^, O = 
(7/A) [exp(-zAZ) - 1], and <I> = (27V A^) [1 - cos(AZ)]. 
In this case the probability of finding this single-photon 



at waveguide site n when launched at /c, can be obtained 
from Pn^k (Z) = (aj^ttn) = \Tk,nf- Equations (g]) clearly 
indicate that the associated probability distribution ex- 
hibits revivals at regular intervals, e.g. at Z = 27ts/X {s 
being integer). At these revival points, for single-photon 
excitation, all the Tk^n coefficients of Eq.(j4]) vanish ex- 
cept for n = /c , i.e., the probability collapses into the 
initial waveguide site k. Figure [T] depicts this process 
for such a Glauber- Fock oscillator array for two different 
values of A. These side views clearly show the previously 
mentioned "collapse" and revivals of the probability at 
Z = 27rs/X irrespective of the site where the photon was 
initially coupled. In all cases, these Bloch-like oscilla- 
tions occur and the photon does not escape into the bulk 
region-towards the right. This in spite of the fact that 
the optical potential linearly increases and the waveguide 
elements get physically closer towards higher values of n. 
We note that unlike standard Bloch oscillations occurring 
in periodic lattices [16-18], in this system the dynamics 
are asymmetric. This broken symmetry is a result of the 
semi-infinite nature of this particular array. 

We now focus our attention on the case where the 
couphng coefficients are Z-dependent, that is when the 
Glauber- Fock oscillator is dynamic. For illustration pur- 
poses we consider the periodic variation: /(Z) = a^o + 
ecos(ti7Z), where /^o is a constant, e is a coupling mod- 
ulation amplitude, and w stands for the modulation fre- 
quency along the propagation direction. In this dynamic 
environment, when a single photon is launched into the 
/c-site, the probability will periodically "collapse" into the 
initial waveguide, at exactly the first zero Z of the func- 
tion |5(Z)p. For the particular example examined here 
these revivals occur when the ratio zu/X = P/Q is a ratio- 
nal number, and P, Q are relatively prime integers. This 
condition is necessary for the two oscillatory processes in- 
dependently occurring in this array to lock together syn- 
chronously. From here one can deduce that Z = 2Pi\ jw. 
This behavior is explicitly illustrated in Fig E] for /^q = 1, 
A = 1, and e = 0.2 for the cases w = 3/4, 2/3, in (a), 
adn (b) respectively. The dashed lines on the other hand 
represent the evolution of the |P(Z)p function which dic- 
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FIG. 1. (Color online) Bloch-like oscillation for a single 
photon propagating through Glauber-Fock oscillator lattices. 
The red arrow indicates the site where the photon is be- 
ing launched. The used lattice parameters are 7=1, and 
A = l/2, A = 4/5 for (a) and (b) respectively. 
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FIG. 2. (Color online) Evolution of the probability distribu- 
tions for single photons propagating through this Glauber- 
Fock oscillator lattices with A = 1 and ii7 = 3/4 for (a), and 
= 2/3 for (b). Dashed lines show the evolution of \B{Z)f 
along Z, predicting the corresponding revival distance at Stt, 
and Gtt, respectively. 



tates the period of oscillations. Note that exact revivals 
do not occur if the ratio tu/A is irrational. On the other 
hand at resonance tu = A, dynamic delocalization occurs. 
In this regime \B{Z)\'^ ex Z, and hence a drift motion is 
induced towards higher site indices [20]. Therefore no 
pure oscillatory behavior is possible at resonance. This 
delocalization process at resonance is depicted in Fig. [H 
where it is shown that the probability of finding the pho- 
ton among the waveguides gradually shifts towards the 
upper side of array oscillator. 

We next consider the quantum dynamics of entangled 
pair of photons, launched either spatially correlated or 
anti-correlated into this class of Glauber-Fock oscillator 
arrays. As we will see, this new class of photonic lat- 
tices can tailor the bunching and antibunching behavior 
of path-entangled biphotons (spatially extended state). 
Conceptually, photon pairs (biphotons) with correlated 
positions will couple into the same waveguide (within a 
certain excitation window W) with an equal probability. 
In this regime, the corresponding input state is written 

2 / , \ 2 / , \ 21 



as 
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|0). 



Throughout our paper, /, and / will represent the first 
and last waveguide site within the excitation window. 
This input state can be generated by placing the waveg- 
uide array immediately after a type I collinear degener- 
ate narrow-band spontaneous parametric down conver- 
sion (SPDC) thin-crystal source[2l|. On the other hand, 
an entangled pair of photons with anti-correlated po- 
sitions corresponds to the physical situation where the 
photon pair is always coupled to waveguides on oppo- 
site sides of the excitation window W-again with equal 
probability. Thus, the input state is written as \^|^A) = 



VVW 



a\al 



|0), where W is 



even and (i?, R') represent the fioor and ceiling integers 
of the quantity (/ + /)/2. In order to obtain the correla- 
tion between the array modes, we analyze at the output 
the coincidence rate at waveguides p and which is given 
by Tp^q = {ajpO^^aqap) . In this case one can show that for 



FIG. 3. (Color online) Theoretical evolution of the probabil- 
ity of having one photon at waveguide n when launched at 
channel 0, Pn,o{Z) — (al^an), under the resonance condition. 



correlated inputs I'^c)^ the correlation map is described 
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der to demonstrate these effects and for comparison pur- 
poses, we will always assume here that = 10 with the 
excitation contained between (/, /) = (0,9). 

The evolution of the correlation map Tp^q as a func- 
tion of distance when f{Z) = 1 and A = 0.5 is then ex- 
amined. In this case, revivals are expected at multiples 
of Z = 47r. When the input state is initially correlated, 
the map fiips and antibunching occurs at Z = tt, and 
Z = Sir while it returns to a broadened bunched state 
in the middle of a cycle. This situation is altered when 
an anti-correlated biphoton input is used. In this case 
the map tends to fiip over to that of a bunched state at 
Z = TT, and Z = 37r while in the middle of this oscilla- 
tion it attains a correlation mixture- with bunching being 
predominant. This evolution is altogether different from 
that occurring in uniform lattices [13j. The present dy- 
namics is a direct outcome of the revivals and of the phase 
acquired upon refiection from the boundary of this semi- 
infinite Glauber-Fock oscillator array-which is absent in 
periodic arrays. We next consider the evolution of corre- 
lations when two periods are simultaneously involved in 
the lattice, e.g. when the function /(Z) is periodic. For 
this example we again take /(Z) = hzo + ecos(tz7Z), with 
A^o = 1, A = 0.5, e = 0.2, and zu = 3/4 in which case the 
revival period is Stt. For a correlated input \^|Jc) the cor- 
relation map exhibits periodic transitions from bunching 
to antibunching FigsHJ^a-e). However, at the half-cycle 
point the bunching is now not entirely complete, FigHJ^e), 
due to the incomplete revival of the single photon trajec- 
tories. This scenario becomes very different when the 
initial biphoton state \iIja) is anti-correlated Fig|4l^f-j). 
The correlation dynamics corresponding to both cases are 
shown up to half a cycle (47r). Evidently, right after the 
origin, bunching is seen to occur (FigHJ^g)) while midway 
in the cycle signatures of antibunching behavior appear. 
This latter pattern is different from that obtained be- 
fore (/(Z) = 1) when only one oscillation frequency was 
involved in the Glauber-Fock oscillator. 
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FIG. 4. (Color online) Left column shows the quantum cor- 
relation maps when the input state lipc) is coupled into a 
modulated-tilted Glauber-Fock oscillator array (A = 1/2, vj = 
3/4). Similarly, right column depicts the correlation evolution 
for the input state I-^a). 



We have also explored the response of this system 
under dynamic delocalization conditions. Figl5] depicts 
again the correlations for the same parameters used in 
the previous case, except that here X = m = 1. In this 
delocalization regime, a correlated input lipc) tends to 
initially oscillate between bunching and antibunching 
as it was shown in Figs. Hl^a-c), and eventually settles 
into antibunched state, FiglSfa). On the other hand, for 
an anti-correlated bi-photon input \^l^A) the entangled 
photons very quickly and irreversibly become bunched 
and they remain in this state, FiglSfb). The reason why 



delocalization itself affects the correlation evolution has 
to do with the fact that in this case the photons tend 
to eventually escape into the bulk of lattice-away from 
the boundary. Simulations indicate that by adjusting 
the two oscillation frequencies one can at will lock the 
output into a particular bunching/antibunching state. 
In essence the presence of revivals of quantum states (or 
absence of revivals) allows one to engineer the quantum 
dynamics in this class of dynamic Glauber-Fock oscilla- 
tor lattices. 
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FIG. 5. (Color online) Correlation matrices at Z = IStt for 
the initial states (a) \ipc), (b) \iPa) and (c) \iPf)- 



Finally, we also consider this same arrangement when 
Fermionic-like input states are used [22]. Figure [5fc) 
shows how the correlation evolves in this case under de- 
localization conditions for parameters identical to those 
used in FiglS] The input in this case is assumed to be of 
the type = b^j^bj |0). For this input, the antibunching 
behavior in the correlation matrix is evident. 

In conclusion we have shown that a new family of dy- 
namic arrays, the so-called Glauber-Fock oscillator lat- 
tices can be used as a way to mold the quantum evolution 
of path-entangled photons. In these systems revivals and 
dynamic delocalization are possible-each leaving a spe- 
cific mark on the correlation map. If the two oscillation 
periods associated with these Bloch-like oscillators are ir- 
rational with respect to each other, the dynamics become 
aperiodic. At this point several intriguing questions re- 
main. For example, of interest will be to examine how 
such structures respond to other maximally entangled 
states (like NOON states) or whether they can be used 
to synthesize other quantum states of interest. The re- 
sponse of these lattices may be also useful in considering 
other classes of problems in other physical configurations 
having similar quantum analogues like those of the Bose- 
Hubbard or Jaynes-Cummings type with time varying 
couplings [Hi [23I . 
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